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THEORETICAL STUDY OF PRODUCTION OF 
UNIQUE GLASSES IN SPACE 

1.0 INTRODUCTION AND SUMMARY 

The overall objective of this program is to develop analyti- 
cal functional relationships describing homogeneous nucleation 
and crystallization in various supercooled liquids. The time 
and temperature dependent relationships of nucleation and cry- 
stallization (intrinsic properties) are being used to relate 
glass forming tendency to extrinsic parameters such as cooling 
rate through computer simulation. Single oxide systems arc being 
studied initially to aid in developing workable kinetic models 
and to indicate the primary materials parameters affecting glass 
formation. The theory and analytical expressions developed for 
simple systems is then extended to complex oxide systems. A 
thorough understanding of nucleation and crystallization kinetics 
of glass forming systems provides a priori knowledge of the ability 
of a given system to form a glass. This will facilitate the dev- 
elopment of improved glasses by providing a firm theoretical/ 
analytical basis for improved manufacturing techniques such as 
in-space manufacture. The ultimate objective of space manufac- 
ture is to produce technically sign 1 ' L 'icant glasses by extending 
the Earth- limited regions of glass formation for certain composi- 
tions, or by achieving glass formation in other compositions that 
are not glass formers based on empirical Earth observations. 

The literature has been reviewed and critically analyzed, 
and kinetic equations have been developed for homogeneous and 
heterogeneous nucleation and subsequent crystal growth. The 
kinetic relationships have been applied to known glass formers 
and non-glass formers. It was found that the models qualitatively 
predict earth-observed behavior for these systems. We can now 
proceed to apply the functional relationships to more complex 
materials with the goal of utilizing the in-space environment 
in producing technically improved glasses. 
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2.0 TECHNICAL DISCUSSION 


The concept of glass forming tendency derives from the 
definition of glass. Morey (1) has defined glass as an inorganic 
substance in a condition which is continuous with, and analogous 
to, the liquid state of that substance; but which, as a result 
of having been cooled from a fused condition, has attained so 
high a degree of viscosity as to be for all practical purposes 
rigid. Hence, glass-forming materials are ones in which there 
is sufficient transient bonding in the melt to produce a highly 
viscous liquid upon cooling. In the general sense, Morey's 
restriction of glasses to only inorganic materials can be re- 
laxed to include organic materials as well. The supercooled 
amorphous state of aggregation of matter comprising a glass is 
unstable relative to the solid crystalline state. Therefore, 
glass forming tendency is related to the mechanisms and para- 
meters that prevent the liquid-solid transformation from occurring. 
Studies of the crystalline transformation can be approached from 
structural, thermodynamic, and kinetic viewpoints. We have chosen 
to adopt a kinetic viewpoint since in general glass formation 
is not related to whether or not a given material can form a 
glass, but rather how fast must the liquid be cooled from its 
fused condition to do so. 

Therefore, we consider a material to be a glass if it can 
be cooled from its liquid state rapidly enough to avoid a certain 
predetermined degree of crystallization. The kinetics of cry- 
stallization of a liquid are determined by two parameters, the 
nucleation rate and the crystal growth rate. The liquid-solid 
transformation occurs by a two-step process of nucleation of 
crystalline embryos and subsequent growth. Nucleation and growth 
rate temperature dependence are illustrated qualitatively in 
Figure 1. Temperature T-^ * T m is the thermodynamic fusion tempera- 
ture where the solid and liquid phases are in co-existence. Above 
T the material is in the liquid phase. As the liquid is super- 
cooled below T , growth can theoretically occur between temperatures 
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1 Figure 1 NUCLEATION AND GROWTH RATE OF CRYSTALS IN GLASS 

AS FUNCTION OF TEMPERATURE 

I 




Room Temp. 


and T^. However, the embryo nucleation that is necessary 
before growth can proceed only occurs between temperatures T2 
and T^. In the temperature region bounded by T2 and nuclea- 
tion and growth occur simultaneously, i.e., conditions are 
favorable for the complete crystallization transformation pro- 
cess, nucleation and growth. Therefore, the (T2 - T^) temperature 
region is most critical to the prediction of glass formation. 

The present program deals with characterizing the nuclea- 
tion and growth characteristics of several materials, and pre- 
sentation in a form similar to Figure 1 . With the analytical 
relationships describing nucleation and crystal growth, computer- 
simulated melt-quench experiments provide information regarding 
the cooling rates necessary for glass formation (i.e., cooling 
rates necessary to pass through the (T2 - T-j) temperature region 
rapidly enough to prevent crystallization). 

The kinetic treatment of glass formation necessarily 
starts with derivations of nucleation and growth expressions. 

In the following sections these parameters are discussed in some 
detail. Succeeding sections will deal with the volume fraction 
transformed and critical cooling rates, initial computer-modeled 
results, and a discussion of the problems in predicting the 
behavior of complex multi- component materials. 
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3.0 KINETICS OF GLASS FORMATION 

Equations describing the kinetics of homogeneous nuclea- 
tion are developed in this section with allowance for transient 
effects. The effect of the presence of nucleating heterogeneities 
are also discussed. Models for subsequent crystal growth are 
presented. Finally, the method of relating the volume frac* . .ji. 
of material transformed to the crystalline state to applies, cool- 
ing rate through computer Simulation is discussed. 

3.1 Theory of Homogeneous Steady State Nucleation 

A supercooled liquid is a metastable phase relative to 
the solid phase as indicated by free energy considerations. The 
liquid system below the fusion temperature Tm tends toward thermo- 
dynamic stability by lowering its free energy through the crystal- 
line transformation. The existence of a supercooled liquid phase 
below Tjj, (i.e., a glass) is the consequence of 1) a surface energy 
barrier between the solid and liquid state, and 2) the kinetically 
inhibited movement of molecules that prevents arrangement in an 
ordered system (i.e., crystalline phase). 

The process initiates as statistical molecular density 
fluctuations causing clustering of molecules or atoms. The 
clusters are called nuclei, embryo particles of solid crystalline 
(transformed) material (transformation to the solid crystalline 
phase can be considered merely as molecular rearrangement into 
an ordered structure). In tending toward thermodynamic stability 
there is a volume free energy decrease as an embrvo is created 
(i.e., in the transformation from liquid to crystalline). However, 
the formation of an embryo means the formation of a boundary, 
the embryo- liquid interface, with a resulting system free energy 
gain due to the interfacial surface energy. Jackson (2) discusses 
the stability of such a nucleus in an undercooled liquid with 
respect to these two phenomena. For small embryos the surface area 
is large, relative to the volume, so that the surface energy domi- 
nates embryo behavior. Small embryo can decrease the total free 
energy of the system (liquid and nuclei) by shrinking to reduce 
their surface area (i.e., dissolving into the liquid melt). 
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Larger size crystalline embryos are controlled by the volume 
free energy term. Large nuclei can reduce the totl free energy 
of the system by growing larger, creating more transformed 
crystalline material (more ordered, more stable). A balance 
between these tendencies defines the critical size nucleus. 

A nucleated embryo smaller than the critical size will dissolve. 

A nucleated embryo larger than the critical size will continue 
to grow. 

Consider a liquid melt at temperature T with order fluctua- 
tions. This system can be described as a steady state concentra- 
tion of ordered regions (crystalline in structure) of various 
sizes. The change in free energy of the system, aF, due to 
a local fluctuation creating a spherical nucleus may be expressed: 

AF " - y - TTr^Af v + 4TTr^Af g (1) 

where r ■ embryo radius 

Af v * free energy difference between the liquid and solid 
state, per unit volume 

Af_ ■ interfacial (surface) free energy between the phases, 
s 

per unit area. 

The system free energy change has a maximum value, aF*, 
for some critical nucleus radius, r*. The critical nucleus 
size, r*, represents the smallest size embryo that can grow with 
I a decrease in free energy to form stable nuclei. The critical 

( nucleus size is derived as follows. Differentiating aF with 

r respect to r gives 

1 

-2AL- - -4nr I 2 Af„ + 8nrAf_ (2) 

f Sr v 8 

Expression (1) has a maximum, aF*, for a value of r ■ r* satis- 
fying 

- o (3) 

Sr 
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Setting the r.h.s. of expression (2) equal to zero and solving 
for r * r* we obtain 


r* - 




(A) 


Evaluating the expression for aF (Equation (1)) for r » r* 
yields aF*, the minimum work required to form a stable nucleus: 


AF 


* 


16tt 

3 



(5) 


The term aF* is the thermodynamic barrier to the nucleation 
process of forming stable nuclei. We now turn our attention to 
obtaining the surface and volume terms in aF* in more useful 
form. 


Turnbull (3) has shown that the liquid-crystal surface 
tension, Afg, can be related to the heat of fusion, Ab^., by the 
expression 


Af 


s 


% n -1/3 v - 2/3 
a m 


( 6 ) 


where Ah f * molar latent heat of fusion of the crystalline phase 

r _ i 

(cal moie ) 

23 - 1 

N ■ Avagadro's number (/^6xl0 molecules mole ) 

V m - molar volume of the crystalline phase (cc/mole) 

Tl dimensionless term a relates the proportionality of Ahf 
and Af c , and is constant for a given type of fluid (3). For 
metals 2, for more complex materials Therefore, 

a has an approximate range 


2< a < 3.3 


(7) 


The parameter a is defined such that physically the recipricol 
of a, 1/a , is equal to the number of monolayers per unit area 
of crystal which would be melted at T m by an enthalpy aH ■ Af g * 

The volume free energy difference between the liquid 
and solid states, A^ , is expressed : 
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-) dT 


which is estimated for small supercooling AT ■ T m - T: 

As f AT 

. r _ -L / r\\ 


where ASj is the molecular entropy of fusion (cal mole'^deg”^) . 

Since the entropy of fusion is related to the heat of fusion. 

Ah, (cal/mole) 

r Ah, 

As, = (10) 


the expression for Af v (P er unit volume) can be expressed 


At large degrees of undercooling AT = - T, this expression 


is modified by the factor T/T 


Af (large AT)= 

V 


T 

( ) 


Limiting our analysis to the case of small undercooling, 
use of Equations (6) and (11) for flf s and ^^respectively, permits 
expression of r* from Equation (4) in more useful terms: 

2T V 

r* - (13) 

aN i/J AT 


For clusters (embryo) below T m with radius r < r*, nuclei form 

and dissolve because the surface energy involved in cluster 

formation is greater than the free energy change accompanying 

solid formation. For r > r*, continuous nuclei growth will oc- 

2 

cur since the surface energy is growing only proportional to r 
while the bulk (volume) free energy term involved in solid for- 

3 

mation is growing proportional to r . Above T crystalline 


1 

I 
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embryo are unstable. Critical cluster radii are illustrated 
in Figure 2 for various oxide systems. 

The thermodynamic barrier to nucleation, AF*, is expressed 
in more useful terms, employing Equation (6) and (11) for Af g and 
Af respectively: 2 

AF* - -MZ ^4 (W) 

J a J NAT 

It will be convenient to express AF* in terms of another dimen- 
sionless parameter, p, defined (9) such that 

As ^ = Nkp = RR (15) 

where k = Boltzmann's constant. 

R = 1.987 cal mole ^deg ^ = gas constant 

The heat of fusion, Ah^, thus becomes 

Ah- = As -T = RpT = NkpT (16) 

r t m T m m 

The parameter p has range 

1 < p < 10 (17) 

with p = 1 for monatomic liquids. More complex structures have 
higher entropies of fusion, with p approaching p * 10. 

The expression for AF* (Equation 14) is expressed in 

terms of p : o 

16 tt kpT 

A F* = — — (18) 

3 a 

Having developed an expression based on thermodynamics 
and free energy considerations for the minimum work required for 
homogeneous nucleation (i.e., the thermodynamic barrier to 
nucleation), we must now derive an expression for the rate of 
homogeneous nucleation. If it is assumed that critical size 
clusters (r = r*) are formed by statistical molecular fluctuations 
the nucleation rate will be proportional to a term involving the 
probability of such a fluctuation, exp (-AF*/kT). The steady 
state nucleation rate involving a Boltzmann- like distribution of 
critical size nuclei is expressed: 
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I o = iW n exp(-£F*/kT) (19) 

where n ■ number of molecules per unit volume in the system 

W 0 0 

(n = — ^ ^ 10^ z ) and is a molecular attempt frequency for 
nucleatTon. 

This result is derived from thermodynamics and thermosta- 
tistics. However, the nucleation process involves molecular 
movement. Therefore, a kinetic or diffusional barrier also 
exists for the nucleation process. The steady state nucleation 
rate, I (nuclei cc'^sec”^), is expressed in general terms 

I Q = iW n exp(-&F*/kT)exp(-AG7kT) (20) 

where AG' is the activation energy for molecular motion (during 
nucleation) across the embryo-matrix interface. Thus the nu- 
cleation rate (steady state) is a function of (1) the number of 
molecular or atomistic units available for nucleation, (2) a 
frequency factor describing how often the molecules attempt to 
jump across the liquid-nucleus boundary, (3) a thermodynamic 
barrier to nucleation, the minimum work required to form a 
stable nucleus, and (4) a kinetic barrier to nucleation, invol- 
ving the activation energy for molecular rearrangement. 

If the molecular motion involved in nucleation is treated 

as an activated process (4), a diffusional rate constant for 

nucleation, D , can be defined such that 
n 

D n = a l V n exp(-AG'/kT) (21) 

where a Q is the molecular jump distance. 

Combining Equations (20) and (21) yields the following 
expression for the steady state nucleation rate: 

nD 

I D = y- exp (-AF*/kT) (22) 
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The problem now arises as how to evaluate the nucleation 
diffusional rate constant, D n . For liquids which crystallize 
without change in composition (single component systems, for 
instance) long range diffusion processes are not required. All 
of the atomic or molecular units required for the ordered cry- 
stalline structure are in the local vicinity of the liquid-crystal 
interface. This type of transformation is termed non- reconstructive. 
The activation energy, AG', for diffusion at the liquid-crystal 
interface will be roughly the same order of magnitude as the 
activation energy for viscous flow. This is the case for a 
liquid transforming to solid without change in composition since 
the movements of the atoms or structural units on the nucleation 
surface is similar to the reorientation of structural units and 
bond switching in the flow of a viscous liquid. Therefore, 
in this simple case, the nucleation diffusional rate constant, 

D n , will be approximately equal to the self diffusion coefficient 
of the undercooled liquid: 

D s < 23 > 

The liquid self diffusion coefficient, D , is related to the 

s 

bulk viscosity, T), through the Stokes-Einstein relation 

D — (24) 

3 rr a o r| 


where a Q is the diameter of the diffusing species. 

Therefore, for the case of a single component substance, 
or a mo.*e complex liquid that crystallizes with no change in 
composition, the steady state nucleation rate is expressed (for 
snu 1 degrees of undercooling below T m ) 


I 


o 


nkT 


exp 


' 16tt P T m 3 

3 

3 ex TAT 


(25) 


where axl terms have been previously defined. 

The general temperature dependence of the homogeneous 
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nucleation rate described by Equation 25 is shown in Figure 1. 

For small degrees of undercooling below T m , AF*, the thermo- 
dynamic barrier, is large since Af v » the volume free energy 
change in the transformation is small. This results in low 
nucleation rates. With further supercooling A^ v increases 
until the kinetic barrier to nucleation. This condi- 

tion results in maximum nucleation rate. For a large aT the 
nucleation rate decreases to a negligible level as AG*> > AF*. 

In the case where a compositional change accompanies the 
liquid- solid transformation, the steady state nucleation rate 
(Equation 25) must be modified to account for the long range 
diffusion processes that are required for molecular rearrange- 
ment. This subject will be treated in Section 6.0. 

3.1.1 Time- Dependent Homogeneous Nucleation 

In this discussion transient effects on nucleation are 
considered. Glass forming tendency is related to how fast 
the system can pass through the region of simultaneous nuclea- 
tion and growth (i.e., the (T^T^) region shown in Figure 1). 

The time spent at any particular temperature level may then be 
less than the time required to establish a steady nucleation 
rate (i.e., the time required to build up the required Boltzmann- 
like distribution of embryos). 

Incorporation of time-dependence into nucleation theory 
has been treated by Hillig (5) and Hammel (4), and takes the 
form 



where = transient nucleation rate 

I = steady state nucleation rate 
t - time 


Hammel (4) has expressed the parameter t (after Collins (6)) 

as: 144kTT_V_ 1 ^ 3 

^7— (27) 


144kTT V 
_ m m 

L “ 7X71 


N a a Q as^aT D 
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where D is the appropriate diffusion coefficient (given by the 

Stokes-Einstein Equation for transformation without change of 

composition) and all other terms have been previously defined. 

This equation was derived for the ideal case of instantaneous 

cooling from T to T. 

m 

Hillig (5) derived an approximate expression for t using 

random walk diffusion theory to describe the mean time "t to build 

a critical nucleus: 9 9 

n V L Z r* z 

t - t ■ 7~T~ (28) 

4DV m z X^ 

m 

where and V m are the mole volumes of the liquid and precipitat- 
ing phase, respectively, D is again the appropriate diffusion 
coefficient (D s D for crystallization without change in composi- 
tion), and X is the mole fraction of the precipitating phase in 
the liquid (i.e., X ■ 1 for crystallization without change in 
composition) . 

Hillig' s t parameter has been shown (4) to more accurately 
describe time- dependence of the nucleation rate, based on cor- 
relation with experimental results. Therefore, we shall employ 
Equations (26) and (28) when simulating nucleation in computer 
generated melt-quench experiments. 

The nature of transient effects on nucleation rate is 
illustrated qualitatively in Figure 3. Due to diffusional effects 
successively longer times are required to achieve steady state 

* 

nucleation at successively lower temperatures. If the hypotheti- 
cal material whose nucleation characteristics are shown in Figure 
3 is cooled from T to T, , detectable nucleation would be avoided 
at any cooling rate.* However, a detectable level of homogeneous 
nucleation occurs between temperatures T^ and In order to 

avoid this nucleation in a melt- quenching experiment the tempera- 
ture region T^ - T 2 must be passed through in less than time t 2 * 

* For many materials T^/^n-^.8 
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If this glass were successfully quenched to temperature Tg, 
for instance, it could be held at this temperature level for a 
period of time tg before detectable nucleation would occur 
(e.g., for annealing). 

As shall be discussed in Section 6.0, the transient 
nucleation formulation provides one possible way of dealing with 
multicomponent systems where long range diffusion processes 
are required for molecular rearrangement. 

3 . 2 Heterogeneous Nucleation 

Thus far we have only treated the case of homogeneous 
nucleation in pure substances free of impurities, insoluble 
particles, etc. Foreign surfaces present in the liquid, such 
as container walls or insoluble particles, however, tend to re- 
duce the barrier to nucleation represented by the surface energy 
between the liquid and solid phases. As a result the critical 
embryo size is reduced and the "supercooling temperature" is 
raised. The "supercooling temperature" is the temperature T 2 
where the first detectable nucleation is observed upon cooling 
from the melt. This effect is shown qualitatively in Figure 4. 

At heterogeneous surface sites the computation of the 
critical embryo size becomes complicated because of the various 
surface energies involved and the particular geometry of the 
embryo surface. Furthermore, the wettability of the foreign 
surface by the liquid must be taken into account. 

For our present purposes we shall only treat the effect 
of heterogeneous nucleation qualitatively, following the treat- 
ment of several authors (references 2,4,5, 7-9) in order to il- 
lustrate the nature of the analytical expressions. 

The work term for heterogeneous nucleation ,AF* 0 , is 
given generally as 

LF* = f (0) . AF* (29) 
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Figure 4 EFFECT OF HETEROGENEOUS NU CLEAT ION ON GLASS FORMING TENDENCY 
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where the factor f(0) is expressed 


f(8) 


[2 + cos9 ) (1 - cos9 V 


The angle 0 is a function of the balance of surface tensions and 
is given by: 

cose - ° sL " 0-8 P (31) 

a LB 


where a gL , ct sp and are the interfacial energies between 
the solid-liquid, solid-impurity surface, and the liquid- impurity 
surface, respectively. 

The lowering of the net interfacial energy in hetero- 
geneous nucleation depends on how well the impurity surface wets 
the nucleating phase in the presence of the liquid. 

Therefore, if heterogeneous nucleation sites are present 
the minimum work required (thermodynamic barrier to nucleation) 
to form a stable nucleus is reduced. The treatment is further 
complicated by the fact t iat the concentration of nuclei de- 
riving from a heterogeneous mechanism will be some function of 
the concentration of impurity sites present, and not given 
simply by a nucleation rate formulation. 

3 . 3 Crystal Growth Kinetics 

Once formed by a homogeneous or heterogeneous nucleation 
mechanism, stable nuclei will continue to grow at a rate determined 
primarily by the rate at which the atoms necessary for growth 
can diffuse to the surface of the crystal, and by the ease with 
which they can free themselves from the attractions of their 
neighbors in the liquid phase and form new bonds in the specific 
positions determined by the structure of the growing crystal (10). 
The fundamental concept here is that growth is considered in 
terms of molecular rearrangement. 

Crystal growth is either interface-controlled or dif- 
fusion-controlled. Interface-controlled growth occurs when 
the rate controlling step occurs at the liquid-crystal interface. 
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This normally occurs for single component glasses, for instance, 
where only short range molecular rearrangement at the liquid- 
crystal interface is necessary. - Diffusion-controlled growth 
occurs when the rate controlling step is the long range diffusion 
of a given species in the bulk liquid. Thus, diffusion-controlled 
growth would normally apply in the case where cry-tallizatirn 
is accompanied by a large change in composition, such as for 
complex multicomponent systems. 

3.3.1 Interface- Controlled Growth Kinetics 

For this work we shall restrict oursel'-e * to the case 
of crystallization without change in composition as we have 
done for the case of homogeneous nucleation in the preceding 
sections. Thus, we will treat in this section only the mechanisms 
and kinetics of interface-controlled growth. 

Two basic mechanisms have been proposed for interface- 
controlled growth: continuous growth and lateral growth (11). 

In continuous growth the crystal interface advances by molecular 
incorporation which can occur with equal probability everywhere 
(except for certain anisotropic effects). Lateral growth occurs 
either by a two-dimensional nucleation mechanism or by a screw 
dislocation spiral growth ramp mechanism depending on the cry- 
stal perfection. 

In systems which crystallize without change in composi- 
tion the nature 'f the liquid-crystal interface strongly influ- 
ences the kinetics and morphology of crystallization. Different 
models for interface controlled growth are each based on a 
different assumption concerning the interface and the nature of 
the sites on the interface where atoms are added or removed. 

Again, it is important to think of crystal growth as merely a 
molecular rearrangement process. 

The general form of the growth velocity for all inter- 
face-controlled growth processes is .* 
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H ■ fa 0 ^g exp(-AG"/RT) [ 1 - exp(-AG/RT)] (32) 

where: f ■ fraction of preferred growth sites on the interface 

(i.e., fraction of sites available for growth, 

0 < f < 1). 

Vg ■ frequency factor for molecular transporc at the 

liquid-crystal interface (during growth). 

a « distance advanced by the interface in a unit kinetic 
o 

process (v^\ molecular diameter). 

£G * free energy change accompanying the l* * .Id* crystal 
transformation 

AG" * free energy of activation or kinetic l er for 
the movement of an atom across the liquid-crystal 
interface during crystallization (i.e., growth). 

The kinetic term aG" involved in growth is not necessari- 
ly equal to (or even the same order of magnitude) as the kinetic 
term involved in nr.cleation, aG‘ (Equation (20)). Growth may be 
governed by atomic movements from a great distance away from 
the liquid-crystal interface. Nucleation (initiation of the 
transformation) is governed by atomic movement relatively close 
to the developing nucleus. 

Assuming that we can treat the molecular movements 
involved in crystal growth as simply activated processes, we 
can define a diffusion coefficient, or rate constant, for the 
molecular rearrangement 

D g “ *o\ exp(-AG"/RT) (33) 

in a manner similar to the case of nvcleation described in 
Section 3.1. The general interface-controlled growth expression 
(Equation (32)) thus becomes: 

U ■ — [l - exp (- aG/RT)J (34) 
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Experience on a variety of pure substances which cry- 
stallize without change in composition has indicated that the 
slow, rate controlling step in the interface growth process is 
the same molecular process that occurs in the bulk liquid above 
the thermodynamic fus-f^r temperature, T m (3). Therefore, in 
this most straight torwarr? of cases, the rate const snt for inter- 
face-controlled growth, Dg, will have approximately the same 
value as the liquid self diffusion coefficient, D , and thus be 
approximated by the Stokes-Einstein relation: 


D * d 
g s 


kT 

3^ 


(35) 


The free energy change accompanying the li< uid-crystal 
transformation, &G, is given by 

LG - &s f AT (36) 


The growth velocity for interface-controlled growth without com 
position change is thus expressed 


u 


fkT “ASfAT 

V [ 1 - exp ( i )] 


3Tra o 11 


(37) 


RT 


or in terms of the 8 parameter (AS£ ■= Nkp * Rp) 

fKT 


u 


3Tta 0 r| 


7~ [ 1 - exp G^£) ] 
T 


(38) 


since 

R - Nk (39) 

We will now briefly treat the various models that have 
been proposed for interface-controlled growth. The models and 
discussion presented have been taken from Uhlmann (12). 

3.3. 1.1 Normal (Continuous) Growth 

In normal or continuous growth atoms can attach to or 
be removed from any site on the interface. Thus, there are no 
preferred growth sites in the interface, and f in Equation (37) 
becomes unity. 
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u 


(40) 


= kT ., 
3tra 0 2 r 1 


[1 


-&s f AT 

exp ( )] 

RT 


For small departures from equilibrium, this model predicts a 
linear relation between growth rate, u, and undercooling, AT. 
For this model to correlate with experimental data the liquid- 
crystal interface must be rough on an atomic scale. 

3. 3. 1.2 Screw Dislocation Growth 


In the screw dislocation model growth occurs at step 
sites provided by screw dislocations intersecting the interface. 
The fraction of preferred growth sites, f in Equation (37), 
is expressed 


f = 


a Q A.h f A.T 

4rrcrT V 
m m 


AT 

2ttT 

m 


(41) 


where it has been assumed that only molecular transport within 

a molecular diameter, a Q , of the dislocation ledge results in 

attachment. Fir small undercooling, &T, this model predicts 

2 

.. growth rate which varies with &T . For the screw dislocation 
model to apply, the interface must be smooth on an atomic scale, 
and be relatively imperfect in the crystallographic sense. 

3.3. 1.3 Surface Nucleation Growth 

According to this model growth takes place at step sites 
provided by two-dimensional nuclei formed on the interface. The 
growth rate can be expressed 

w r exp (-B/TAT) (42) 

where the exponential constant is 

Tra V T Op 2 

B ° m m E (43) 

kAhf 


The term Og is the specif i edge rurface energy of the nucleus. 
The frequency factor, can be derived from the growth rate 
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constant : ^ 

§- (44) 

a 

o 

The growth velocity predicted by this model should vary exponenti- 
ally with the undercooling, &T, and for small AT should be unob- 
servably low. For this model to correspond to experimental ob- 
servation the interface must be smooth on an atomic scale and 
be defect-free. 

This comment and the foregoing discussions regarding 
expected morphology for each of interface-controlled growth 
models will have application in the experimental stages of our 
NASA program (to be discussed later). For materials which have 
not been extensively studied, morphological studies (SEM, etc.) 
will indicate, by inference, which growth model applies to a 
given material. Utilizing empirical results of melt-quench 
experiments to correlate with our kinetic relationships will 
help to indicate how space processing can be utilized to pro- 
duce glasses unobtainable on Earth. 

3.3.2 Significance of Entropy of Fusion 

As discussed by Uhlmann (12) growth rates depend critical- 
ly on the molecular structure of the liquid-crystal interface. 

The interface structure, in turn, depends significantly on a 
bulk thermodynamic property, the molecular entropy of fusion, 

A s ^ . 

For materials characterized by low entropies of fusion 
(ASj < 2R) the liquid-crystal interface should be rough on an 
atomic scale, defects should not affect growth, and normal growth 
kinetics are predicted (i.e., Equation (40)). On the scale of 
light microscopy the crystal- liquid interface should be non- 
faceted. 

For materials characterized by large entropies of fusion 
(As f > 4R) faceted interface morphologies should be observed, 
defects should be important to growth, and the kinetics of the 
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normal growth model should not apply. Except as limiting cases, 
the observed growth velocities for large As^ materials should 
not agree well with behavior predicted by screw dislocation 
model or the surface nucleation model either. 
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4.0 VOLUME FRACTION TRANSFORMED AND CRITICAL COOLING RATES 


Glass forming tendency can quantitatively be describe! 
in terms of the volume fraction of crystalline material formed 
during a certain quenching operation. Uhlmann (13) has expres- 
sed the volume fraction, V£, crystallized in time, t, (for small 
Vf) as 

v f ^1/3 ttI u 3 t 4 (45) 

where I = nucleation rate (nuclei cc’^sec”^) 

U * growth velocity (cm sec”'*') 
t s time (sec) 


This relation is valid for interface-controlled growth, and 
thus applies to single component or congruently melting (i.e., 
without change in composition) compounds. 


The critical cooling rate for glass formation, i.e., the 
minimum cooling rate necessary to avoid a certain volume frac- 
tion of transformed material, can be estimated by a procedure 
employed by Uhlmann (13). T-T-T (time -temperature- transformation) 
curves are constructed as illustrated qualitatively in Figure 5. 
Using Equation (45), the time required to transform a given 
amount of material is calculated as a function of temperature. 

A volume fraction transformed v^ * 10”^ is regarded as just- 
detectable. The data are presented graphically as AT vs t in 
the T-T-T plot. The critical cooling rate can be approximated 
(13) by the expression 


dT 

“cTt~ 


critical 


^N 

T N 


(46) 


where AT N ~ T m - T N 

T^ ■ temperature at the nose of the T-T-T curve 
= time at the nose of the T-T-T curve 
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Figure 5 TIME-TEMPERATURE-TRANSFORMATION CURVE 
(TAKEN FROM UHLMANN (13)) 
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5.0 COMPUTER MODELED RESULTS FOR SELECTED SYSTEMS 

The kinetics developed for nucleation and growth are 
now applied to several materials systems. The well-characterized 
silica, Si 02 , system is investigated for the purpose of deter- 
mining the relative importance of the various kinetic and thermo- 
dynamic parameters in controlling transformation behavior. The 
I^O^ system is studied to determine if the derived kinetic re- 
lationships will predict its good glass forming qualities, and 
imply a physical reason for the observed behavior. The alumina, 
^*2^3) s y stem studied to determine if the kinetics can pre- • 
diet the poor glass forming tendency that is observed on earth, 
and imply a physical reason for this phenomenon. 

5.1 The Silica System (SiC^) 

The pertinent materials properties for the silica system 
to be used in the derived kinetic nucleation and growth expres- 
sions are tabulated in Table X. Viscosity data for the silica 
system are shown in Figure 6, compared with viscosity data for 
other oxide systems. 

These materials parameters were inserted into the derived 
nucleation and growth kinetics (Equations (25), (26), (28), 
and (37^), with results illustrated in Figures 7 and 8. The 
low temperature cut-off in the nucleation rate has been discus- 
sed previously. Figure 7 indicates the relative time scale of this 
effect. A nucleation rate of 1 nuclei cc ^sec ^ has arbitrarily 
(but by popular custom (5)) been taken as the level of detectable 
nucleation. As illustrated in Figure 7, detectable homogeneous 
nucleation for SiC^ does not occur upon cooling from the melt 
(T m ) until roughly 1700° to 1750°K. Above about 1200°K the 
thermodynamic barrier to nucleation controls behavior. Below ^ 

* 1200°K the kinetic (diffusional) factor dominates. The peak 

* growth rate for SiC^ (Figure 8) is roughly 40 A/min, and occurs 

j. at roughly 50 °K undercooling below T ra . The crystal growth rate 

is seen to be zero at T m » the thermodynamic fusion temperature. 

J For small values of supercooling (temperatures during quench 
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TABLE I 

MATERIALS PARAMETERS FOR SlO, 


Ahf = 2000 cal mole" 

T = 2000°K 

m A U 

As^ = ip = 1 cal mole K 

m 

ASr 

ft - -/~.5 

o 

a Q = 2.5 A (vhO-O distance in SiO^ tetrahedron) 

f =1 (normal growth assumed) 

X * 1 (pure substance) 

V - V T = 27.6 cc mole' 1 

m L 

T| (viscosity) : see Figure 6 

a =2.5 (representative value) 

23 -1 

N = 6x10 molecules mole 

n = ~ — = 2x10^ molecules cc”^ 

v m 

AT = T m -T 

k = Boltzmann's constant = 1. 38xl0~^erg c”^ = 3.3xl0"^^cal c”^ 

R = gas constant * 1.987 cal mole - ^deg - ^ 
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Figure 8 CRYSTAL GROWTH RATE (SiO,) 
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below T m ), the free energy term is dominant and increases with 
increasing supercooling (i.e., the growth rate increases with 
decreasing temperature). At large degrees of supercooling 
below T m , the kinetic factor (diffusion) begins to control, 
and the growth rate begins to decrease with decreasing tempera- 
ture. At temperatures far below T m the diffusion kinetics 
dominate: the extremely high viscosity inhibits crystal growth, and 

the growth rate drops to zero. 

The critical temperature region for glass forming ten- 
dency, the region of simultaneous nucleation and growth, is 
illustrated qualitatively in Figure 9 in a "Tg-Tj 1 ' region plot. 

If homogeneous nucleation only is considered, it is observed 
that Si (>2 glass formation depends on the ability to pass through 
the temperature range 1600-1700°K rapidly. The transient 
behavior shown in Figure 7 indicates that this critical region 
must be traversed in a time less than approximately 5x10 sec 
to avoid the liquid-crystal phase transformation. This time 
scale limit is readily attained by commercial earth processing 
methods . 

5 . 2 Hypothetical Glass- Parametric Study 

The computer software that has been developed to generate 
nucleation and growth behavior provides a convenient means of 
determining which materials parameters are of primary importance 
in determining glass forming tendency. This information is 
critical to our overall task of investigating the glass forming 
tendency of unique systems, or in synthesizing systems for space 
manufacture . 

For this purpose a hypothetical glass was postulated 
that possessed the viscosity-temperature behavior of SiC^, and 
the parameters a and |3 were varied over realistic ranges (i.e., 
2^a;<3, 1 < b < 10). The resulting parametric graphical 
representations are presented in Figures 10 and 11, and serve 
to indicate the importance of these parameters on the (steady 
state) nucleation rate. The horizontal line is the level of 
detectability of nucleation (assumed to be 1 nucleus per cm per 
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Figure 9 T^-T^ REGION NU CLEAT ION AND GROWTH PLOT FOR SILICA 
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Figure 10 STEADY STATE NU CLEAT ION RATE FOR A HYPOTHETICAL GLASS (T m - 2000°K 
P « 1, a - VARIABLE); THERMODYNAMIC AND KINETIC BARRIERS 
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Figure 11 STEADY STATE NUCLEATION RATE FOR A HYPOTHETICAL GLASS (T m - 2000°K 
a = 2, B - VARIABLE); THERMODYNAMIC AND KINETIC BARRIERS 
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sec). The temperature at which the nucleation curves intersect 
the horizontal level of detectability line is temperature T 2 
shown in Figure 1. Above T 2 no nucleation should be observable. 
Since nucleation and growth occur simultaneously between T 2 and 
Tp we are interested in decreasing (T 2 -T^); therefore, decreasing 
T 2 will enhance glass forming tendencies. Additionally, we are 
interested in material parameters that tend to decrease the 
slope of the nucleation rate- temperature relationship in the 
region just below temperature T 2 (this is important due to I imita- 
tions in cooling rates attainable). 

Figure 10 illustrates the effect of a for the hypothetical 
glass with the S parameter fixed at unity. It is observed that 
decreasing a decreases T 2 , and thus tends to increase glass 
forming tendency. 

Figure 11 illustrates the effect of p for the hypothetical 
glass, with the a parameter fixed at a value of a = 2. It is 
observed that the g parameter has a larger effect on the nuclea- 
tion rate than does the a parameter. Temperature T 2 decreases 
with increasing p. 

It is apparent that for enhancing glass forming ten- 
dencies large p and small a parameters are required. This is 
based on the observed decrease in the temperature T 2 and de- 
crease in the slope of the nucleation curve near T 2 « This can 
also be seen by considering the expression for AF*, the free 
energy of a critical embryo, which represents the magnitude of 
the barrier to nucleation. A large p represents a material with 
a large heat of fusion; a complicated structure with a large 
amount of energy involved in the phase change. A small a 
represents a material with a large heat of fusion relative to 
the solid- liquid interfacial free energy. Values for p (or ASf) 
for several systems are shown in Table II. 

Thus far we have only discussed thermodynamic parameters 
a and 0. Since ther also exists a kinetic barrier to nucleation 
and growth it would be expected the viscostiy, T|, would be a 
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TABLE II 


PERTINENT 

PROPERTIES OF 

VARIOUS MATERIALS 

T ( rt K) 

Material 

As f 

(cal/mole°K) 

Ah f 

Ical/molel 

0 

Si0 2 

1 

2000 

.5 

2000 

B 2°3 

8.2 

5900 

4.1 

723 

Ge(>2 





a-Al 2 0 3 

11.3 

26,000 

5.6 

2300 

lithium silicates 

5 




potassium silicates 

12 




Ca0.Al 2 0 3 .2 Si0 2 

16.1 




Soda-Lime glass 

8 




Metals 

2.3 
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parameter of primary importance. Uhlmann (13) has discussed 
the importance of two kinetic parameters on the tendency of a 
material to be a glass former: ( 1 ) a high viscosity at the 
fusion temperature, T , and (2) a rapidly rising viscosity with 
decreasing temperature below T m (i.e., large dTi/dT). The latter 
is related to the position of T m relative to the glass transition 
temperature, T . We will not elaborate on the parameters “n and 
d'H/dT in a parametric study. However, they will be considered 
in subsequent discussions of the glass forming tendencies of 
B 2 O 3 and A^O^. 

5.3 The B 0 O 3 System 

In this section nucleation and growth kinetics are ap- 
plied to the B 2 O 2 system to predict the observed good glass 
forming quality for this material, and to imply a physical reason 
for this phenomena. Qualitative indications of glass forming 
tendency can be obtained by considering only steady state nuclea- 
tion, ignoring for the moment transient effects. Steady state 
data represent an upper bound for nucleation behavior. 

The pertinent material properties for the B 2 O 3 system 
are tabulated in Table III, The viscosity of B 2 O 2 is illustrated 
in Figure 12. Over the temperature range 250-500°C the litera- 
ture data shown were averaged and fit to a 6 th degree polynomial 
expression: 

log T| = AT 6 +BT 5 +CT 4 +DT 3 +ET 2 +FT+G (47) 

with resulting coefficients: 

A = 2.24xl0" 14 
B =-4.56x10 
C = 3.82xl0“ 8 
D =-l. 73xl0~ 5 
E = 4. 81xl0 -3 
F - -8. 9xl0 _1 
G = 9.77X10 1 
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TABLE III 

MATERIALS PARAMETERS FOR B,^ 

Ahf = 5900 cal mole”*" (reference 21) 

T m * 450°C * 723°K 

Ahf _ I _ -i 

ASf = y — = 8.2 cal mole K 



a Q “ 2.5 A (assumed) 

f =1 (assumed normal growth) 

X =1 (pure substance) 

a: variable, 2 < a < 3 
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J The computer-generated viscosity data for are a ^ s0 s ^ own 

in Figure 12 for comparison. 

j These material properties were inserted into our computer 

software and the steady state nucleation rate for computed 

1 and plotted for variable a as illustrated in Figure 13. The 

computed growth rate for iH ustrate ^ i n Figure 14. Glass 

forming tendency for is deduced from considerations of the 

' nucleation and growth behavior shown in these two Figures. 

Referring to the growth curve (Figure 14), it is observed 
that the growth rate for is zero at T m , peaks at roughly 700° K, 

and decays to zero by roughly 625°K. The steady state nucleation 
rate for (Figure 13) was computed for the known P parameter 

(R = 4.1) and for a probable range of the a parameter. It is 
observed that for all values of a, the nucleation rate is always 
below the detectable limit (assumed 1 nucleus per cc per sec), and 
thus effectively zero. Therefore, although exhibits a theo- 

retical growth rate much higher than SiOg (see Figure 8), critical 
j size crystal embryo are never homogeneously nucleated. These 

analytical results predict , therefore, the would be an excel- 

l lent glass former since it can not be made to nucleate and crystal- 

■ lize when cooling from the melt, regardless of the cooling rate. 

This prediction is supported by empirical evidence indicating that 
crystallization of from a dry melt has never been observed. (12) 

The parameters of primary importance in describing glass 
forming tendency from a kinetic standpoint are 1) the magnitude 
of the viscosity at T m , 2) the slope of the viscosity-temperature 
function below T m , dri/dT, and 3) the magnitude of the molecular 
entropy of fusion Asf (related to AHf and P). The former two are 
related to the diffusional barrier and the latter to the thermo- 
dynamic barrier to the phase change process, liquid state to 
crystalline state. The magnitude of the melting point viscosity 
for ?>2®2 two orc * ers magnitude less than that of SiC >2 (see 
Figure 6). This would indicate that nucleation-crystallization 
should occur more readily for (lower diffusional, kinetic 
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Figure 14 CRYSTAL GROWTH RATE FOR B,0. 
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! barrier) . Since this is not the case, the factors inhibiting 

I the phase change in are the high molecular entropy of 

j fusion ( 8.2 for roughly 1,0 for Si 02 ), and the steepness 

of the viscosity- temperature relation below T m for (dri/dT 

greater for than for Si 02 as shown in Figure 6 ). 

To illustrate the relative importance of £s£ and dn/dT 
in inhibiting nucleation and crystallization in the s y stem > 

we have used our computer system to describe the nucleation 
behavior of a hypothetical liquid with the viscosity- temperature 
relation of I^Oj, but with the molecular entropy of SiC^ (i.e„, 

ASf = 1), Our hypothetical glass then exhibits the kinetics of 
a material that will not nucleate and crystallize (I^O^) and 
the thermodynamics of a material that will nucleate and cry- 
stallize (SiC^), The nucleation behavior of this hypothetical 
material is shown in Figure 15, with the behavior of 
(a,S£ “ 8,2) also shown for reference. It is observed that the 
nucleation rate of the hypothetical liquid is now above the 
detectable limit. Thus by lowering the thermodynamic barrier 
to the phase change for (i.e,, &S£ or p) we are able to 

produce a hypothetical glass that Wj. 11 nucleate. The growth 
rate for this hypothetical liquid was computed and is shown in 
Figure 16, Comparing the relative positions of the nucleation 
and growth curves for this hypothetical liquid (Figures 15 ar.d 
16) indicates its glass forming qualities. For convenience 
these nucleation and growth curves are qualitatively sketched 
in a "T 2 -T 3 region" plot in Figure 17 It is observed that a 
small region of simultaneous nucleation and growth exists between 
625 and 650 °K„ Thus crystallized could be obtained upon 

cooling from the melL (by holding the liquid between 625 and 650 °K) 
if the molecular entropy of fusion (or heat of fusion) could 
be reduced substantially. This hypothetical example is pre- 
sented as an indication of how we might eventually be able to 
synthesize materials to exhibit desired nude? tion-crystallization 
tendencies , 
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5.4 The AI 2 O 3 System 

Nucleation and growth kinetics also have been applied 
to the A^O^ system. This system exhibits very poor (non-existent) 
glass forming tendencies by conventional earth methods. The 
aim here is to predict this poor glass forming tendency of A^O^ 
and imply a physical reason for this phenomenon. 

The pertinent materials properties for the alumina 
system are tabulated in Table IV. The major difficulty in des- 
cribing nucleation and growth kinetics for the A^O-^ system is 
the complete lack of viscosity data below the fusion temperature, 

T . This situation exists since Al 0 Oo has never been observed 
sn earth as a glass (i.e., at temperatures below T m ). Above 
the melting point, however, data have been reported, ^ a nd are 
illustrated in Figure 6 . With no reference pcint such as the 
glass transition, T^, available, the problem remains as to the 
shape of the viscosity-temperature curve for A^O^ below T m « 

For our purposes we have chosen two cases: 1) extrapolating 
below T m with the shape of the B 2 O 3 curve (i.e., steep dri/dT) 
and 2 ) extrapolating below T m with the shape of the SiC^ curve. 
These two cases thus represent the upper and lower bounds for 
dri/dT of the various oxides shown in Figure 6 . For both cases 
the extrapolated curves were expressed algebraically using curve 
fitting techniques for use in our computer software. 

The homogeneous nucleation and growth characteristics 
of AI 2 O-J are summarized in the "Tj-T^ region" plot shown in 
Figure 18. Behavior depicted by solid lines represents the 
case where A^O^ is assumed to have a dri/dT behavior below T^ 
similar to that of SiC >2 (see Figure 6 ). The dotted lines in 
Figure 18 represent the case where the viscosity-temperature 
dependence of is assumed below the A^O^ fusion temperature, 

T . Since Al o 0 Q has such a low viscosity at its fusion tempera- 
ture (melting point), theoretical growth rates are much higher 
than for either pure SiC^ or Both of the AlpO^ growth 
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TABLE IV 


MATERIALS PARAMETERS FOR AloOo 


T - 2300°K 

m 

o 

a Q ■ 2.5 A (assumed) 
f -I (assumed) 

Ahf ■ 26,000 cal mole”^ (reference 21) 
Ahf _ i _ i 

As^ • ^ — * 13 11.3 ceI mole K 
m 

As f 

(3 " g" “ 3.65 

a - 2, 2.5, 3 (variable) 
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CHARACTERISTICS 


rate curves shown in Figure 18 have peak magnitudes of w^lO^® 

O O / o 

A/min, compared to 40 A/min for Si02 (Figure 8) and 2x10^ A/min 
for 1^3 (figure 14). To illustrate the very large theoretical 
growth rate of the peak growth of (2x10^ A/min) is 

reached by the alumina system within v* 10"^ °K undercooling below 
T m . The 40 A/min peak growth rate of Si02 is reached by A^O-j 
within ^10 ® 6 K undercooling below T m . 

At large degrees of undercooling below T , in the dif- 
fusion-controlled portion of the growth curve, the steepness of 
the viscosity-temperature relation is seen to determine the tem- 
perature at which the growth rate drops to zero. Since we are 
dealing with the approximate upper and lower boundaries of the 
possible viscosity relation of alumina below T , this temperature 
(T^ in a conventional '^-T^ region" nucleation and growth plot) 
lies between approximately 1300° and 2050°K. 

The steady state homogeneous nucleation rate for A^O^ 
is also shown in Figure 18 (for a variable). For the condition 
where the A^O^ viscosity has the general temperature dependence 
of ^ 2^3 below T m (steep cFl/dT) it is observed that homogeneous 
nucleation is below the detectable limit. Where A^O^ is assumed 
to behave more like SiC^, which is more probable, it is observed 
that a detectable level of homogeneous nucleation occurs for <x=3 . 
The shaded region shown in Figure 18 thus represents the region 
of simultaneous (homogeneous) nucleation and crystal growth for 
A^O^ system. 

If homogeneous nucleation and subsequent growth are 
considered, we have shown that to avoid crystallization upon 
, cooling from the melt the temperature region 1300° to 1625 °K 

j must be passed through rapidly. In practice, however, A^Og 

I 

j is known to crystallize almost immediately upon cooling below 

■ the fusion temperature, T m (2300°K). The reason for this dis- 

1 crepancy is that our analytical predictions of glass forming 

tendency have not accounted for heterogeneous nucleation. In 
the case of heterogeneous nucleation any insoluble inpurity or 
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external surface will serve to lower the size of a critical 
embryo and thus effectively increase temperature T 2 (the tempera- 
ture where detectable nucleation first appears upon cooling from 
the melt). This effect is shown in Figure 19. Temperature T 2 , 
corresponding to homogeneous nucleation, is shifted to a higher 
temperature ( 12 ') if heterogeneous nucleation is possible (i.e., 
if extrinsic nucleation sites are present). As illustrated in 
Figure 19, simultaneous heterogeneous nucleation and growth 
can (qualitatively) occur at temperatures just below the fusion 
temperature, T . Most investigators of crystallization phenomena 
indicate that nucleation appears to initiate heterogeneously 
in most materials. Thus an extrinsic property controls the glass 
forming tendency of A^Og, at least on earth where the. complete 
elimination of heterogeneous nucleation sites i$ not possible. 
Complete elimination of external nucleation sites in AlgOg 
is necessary if A^Og glass is to be obtained since the growth 
rates are so high. Therefore, A^Og might be a good candidate 
for space manufacture since processing could be performed con- 
tainerless, with no external surfaces acting as nucleation sites. 
Before the space-processing candidacy of A^Og is determined, 
however, several areas must be investigated in more detail. 

These include: 1) glass quality regarding crystal size and 
concentration if homogeneous nucleation only is possible (i.e., 
considering the high growth rates and attainable quench rates 
from 1625° to 1300°K, perhaps enough crystalline phase would 
be nucleated homogeneously to yield a poor quality glass even 
with space processing), and 2) insoluble impurity levels attain- 
able in Al 2 0g precursor materials (i.e., perhaps enough impurity 
sites will be available for heterogeneous nucleation to make 
elimination of surface nucleation sites (crucible wall) through 
space processing only a second order improvement). Clearly, 
if space processing is to be employed to eliminate heterogeneous 
nucleation sites leading to high quality glasses, then the mechanisms 
of heterogeneous nucleation must be well known before candidate 
materials can be chosen with any confidence. 
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Figure 19 EFFECT OF HETEROGENE 01 ' NUCLEATION ON GLASS FORMING TENDENCY OF Al 2 0 


6.0 APPLICATIONS TO MORE COMPLEX SYSTEMS 


The nucleation and growth kinetics that have been devel- 
oped in Section 3 apply in general to single component substances 
or congruently melting (i.e., without composition change) com- 
pounds. The liquid-solid transformation process for these materi- 
als is termed non- reconstructive. No interatomic bonds within 
the participating molecules need be broken, and only short range 
diffusion processes are required for the transformation to occur 
at the liquid-crystal interface. In this instance, the mole- 
cular movements are treated as simply activated processes, with 
rate constants approximately equal to the coefficient of self- 
diffusion in the bulk liquid. This is convenient since the self- 
diffusion coefficient is related to the viscosity through the 
Stokes-Einstein equation, and viscosity data are more readily 
attainable kinetic data than diffusion data. 


In multicomponent systems, however, the liquid to solid 
transformation often involves bond breaking and/or long-range 
diffusion processes, as discussed previously. In this instance 
the transformation is termed reconstructive. In network liquids, 
such as SiC^, interatomic bonds must be broken in the network 
prior to molecular rearrangement. Since this bond breaking in 
network liquids must also preceed viscous flow or self-diffusion, 
the free energy of activation for these processes will also be 
applicable as the kinetic barrier to the phase charge process (25). 
Thus, in the reconstructive transformation of a network liquid 
we can approximate the activation energy and the rate constant 
in the same manner as for non- reconstructive crystallization: 


AG' (nucleation) = AG" (growth) = AG (flow activation energy) 

CL 


and, 


D * D = D e 
g s 


kT 

3TTa h 
o 


(49) 


(48) 


However, in reconstructive crystallization where a large 
change in composition is involved, long range diffusion processes 
are required to bring the appropriate atomic species to the liquid- 

1 1 T RESEARCH INSTITUTE 


54 



crystal interface. In this case the rate limiting step may be 
the diffusion of a particular (i.e., the slowest moving) species 
in the liquid matrix. In dealing with nucleation and growth 
kinetics and subsequent glass forming tendency of such materials, 
several methods have been postulated. Uhlmann and Chalmers (8) 
suggest that for nucleation involving large changes in composi- 
tion that the kinetic barrier to nucleation, aG' in Equation (20), 
should be taken as the activation energy for diffusion of the 
slowest moving component in the matrix, and that the pre-exponential 
factor, nV , should be reduced by the mole fraction of the pre- 
cipitating component. This presents a problem since the appropri- 
ate diffusion data are not as readily available as bulk liquid 
viscosity data. Furthermore, the kinetic term for nucleation, 

&G', is not necessarily equal to or even the same order of mag- 
nitude as the kinetic term for growth, aG" in Equation (32). 

This is due to the fact growth is governed by atomic movements 
from a great distance from the interface, whereas the molecular 
movement involved in the nucleation process occurs relatively 
near the liquid-crystal interface. In general, growth in multi- 
component systems where composition is a variable tends to be 
diffusion-controlled rather than interface-controlled. Hammel (4) 
has discussed the formulation for computing the volume fraction 
of material transformed for diffusion controlled growth: 


_ 8tt o3 n ,3/2 T .5/2 
~U~ S D I o t 


(50) 


where D' is the diffusion rate constant and S is a supersatura- 
tion term. Again, however, diffusion coefficients are required 
that are not readily available. 


The question arises, then, of how we can predict the 
glass forming tendency in systems where diffusion controlled 
behavior is expected due to large compositional changes during 
the transformation, and for which the required diffusion data 
are not available. Hillig (11,26) has proposed that incorpora- 
tion of a transient nucleation rate (5) will adequately account 
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for long range diffusional effects that occur in a reconstructive 
transformation. Hammel (4) has applied transient nucleation 
analysis (Equations 26, 27, 28) and diffusion controlled growth 
analysis (Equation 50) to predict the volume fraction of cristo- 
balite (Si02) precipitating from E glass (13Na02* llCaO* 76Si02) 
during cooling from the melt. In the nucleation rate expression 
Hammel used a value of X = .1 (mole fraction of precipitating 
phase in the melt. Equation 28) to account for the long range 
molecular rearrangement. A value of X equal to unity is employed 
for a pure single component substance or a congruently melting 
compound. The determination of the appropriate value for X to 
use in a given reconstructive transformation is not clearly 
defined, however. Hammel (27) and Hillig (26) have proposed 
that X be determined by considerations of 1) what the expected 
rate- limiting species will be and 2) its concentration in the 
melt. In Hammel' s treatment of diffusion controlled growth 
described above, a value of S in Equation (50) was determined 
using Frank's (28) tabulated values for diffusion controlled 
growth. 

Unfortunately, the state-of-the-art in this nucleation- 
crystallization area has not reached the point where predicted 
transformation kinetics correlate well with reality in all in- 
stances. Especially in the case where large composition changes 
accompany the transformation, we may have to resort to empirically 
derived nucleation and growth kinetics as discussed by Uhlmann (12). 
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7.0 CONCLUSION AND FUTURE WORK 

The emphasis in this program is in developing the analyti- 
cal tools to permit prediction of glass forming tendency in uni- 
que oxide systems. Unique systems, for our purposes, means 
systems that are not good glass formers by conventional Earth- 
processing means. 

A fairly comprehensive treatment of nucleation and growth 
kinetics in pure substances has been presented. The derived 
transformation kinetics have been successfully applied to a well- 
characterized system (Si 02 ), an excellent glass former 
and a poor glass former by conventional means (A^O^). The 
kinetic and thermodynamic parameters of viscosity and entropy 
of fusion have been shown to be the primary materials parameters 
controlling glass forming tendency. 

For complex multicomponent systems where diffusional 
effects predominate, the state-of-the-art is not nearly as far 
advanced as for simple substances. The transformation kinetics 
of materials which crystallize with a large compositional change 
are most probably governed by the long range diffusion of the 
slowest moving species. The general lack of specific diffusion 
data, however, dictate that simplifying assumptions be made, 
such as the validity of the Stokes-Einstein equation relating 
bulk diffusion and viscosity. 

With the nucleation and growth kinetics that have been 
developed, we can now begin to determine how a processing tech- 
nique such as space processing can be utilized to produce tech- 
nically significant glasses unobtainable on Earth. Our work 
has lead us to the point where we need to achieve a better under- 
standing of real, nonsimple systems. This mainly involves a 
better understanding of diffusion-controlled kinetics and effects 
of heterogeneous nucleation phenomena. Since adequate literature 
data to be used in our analytical relationships are non-existent 
we must rely on empiricism. A series of actual melt-quench 
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experiments in the laboratory will permit us to correlate the 
behavior of potential space processing candidates with the 
theory ard analytics that have been developed. It is believed 
that the coupling of our kinetic equations with empirical evi- 
dence and post-quench examination will permit reasonably accurate 
predictions of the glass forming tendency in unique systems. 

Our goal is to be able to predict the level of product improve- 
ment that is obtainable through in-space materials processing. 

To achieve this objective we shall employ our idealized 
kinetics to try to extend the glass forming region of calcium 
aluminate, Y2O3, 0320^, and mullite, for instance. The subject 
of glassy mullite compositions is currently receiving attention 
in the literature, and calcium aluminate has application (for 
high alumina concentrations) in ir-transmitting systems (i.e., 
sapphire is a well known ir-transmitting material.). 

To facilitate the correlation of theory and experiment 
we are currently working with the North American Rockwell group* 
whose program entails the production of unique glasses by a laser 
spin melting technique. In this cooperative effort we shall 
conduct a series of experiments designed to aid in the determina- 
tion of candidate space processing materials, and to predict 
the level of product improvement to be expected in space manufac- 
ture. Soecifically, the experimental observations will help us 
to fit heterogeneous nucleation and diffusion controlled behavior 
into our kinetic equations in the absence of the required litera- 
ture data (such as surface energies, wettability, and diffusion 
coefficients). Areas of investigation in this cooperative effort 
include 1) degree of superheat (above T m ) necessary to eliminate 
volume heterogeneous nucleation (impurity) sites, 2) effectiveness 
of containerless processing in eliminating surface heterogeneous 
nucleation sites, 3) cleanliness of the sample environment (at- 
mosphere) required in space processing of unique materials, etc. 

* NASA Contract NAS8-28991, R. A. Happe, Principle Investigator. 
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The results of this work will appear in subsequent progress 
reports . 
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